AFRL-RY-WP-TP-2008-1160 


COMPARATIVE  STUDY  OF  COHERENT,  NON¬ 
COHERENT,  AND  SEMI-COHERENT  INTEGRATION 
SCHEMES  FOR  GNSS  RECEIVERS  (PREPRINT) 

Chun  Yang,  Mikel  Miller,  Erik  Blasch,  and  Thao  Nguyen 
Sigtem  Technology,  Inc. 


APRIL  2007 


Approved  for  public  release;  distribution  unlimited. 

See  additional  restrictions  described  on  inside  pages 


STINFO  COPY 


AIR  FORCE  RESEARCH  LABORATORY 
SENSORS  DIRECTORATE 

WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7320 
AIR  FORCE  MATERIEL  COMMAND 
UNITED  STATES  AIR  FORCE 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

2.  REPORT  TYPE 

Conference  Paper  Preprint 

4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 

COMPARATIVE  STUDY  OF  COHERENT,  NON-COHERENT,  AND  SEMI-  FA8650-05-C- 1 828 

COHERENT  INTEGRATION  SCHEMES  FOR  GNSS  RECEIVERS  (PREPRINT)  5b.  GRANT  number 

5c.  PROGRAM  ELEMENT 
NUMBER 

65502F _ 

5d.  PROJECT  NUMBER 

3005 _ 

5e.  TASK  NUMBER 

13 _ 

5f.  WORK  UNIT  NUMBER 

30051 3  CY _ 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION 

Reference  Systems  Branch  (AFRL/RYRN)  REPORT  NUMBER 

RF  Sensor  Technology  Division 
Assessment  and  Integration  Branch  (AFRL/RYAA) 

Sensor  ATR  Technology  Division 
Air  Force  Research  Laboratory,  Sensors  Directorate 
Wright-Patterson  Air  Force  Base,  OH  45433-7320 
Air  Force  Materiel  Command,  United  States  Air  Force 
9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 
Sensors  Directorate 

Wright-Patterson  Air  Force  Base,  OH  45433-7320 
Air  Force  Materiel  Command 
United  States  Air  Force 


10.  SPONSORING/MONITORING 
AGENCY  ACRONYM(S) 

AFRL/RYRN 

11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER(S) 

AFRL-RY -WP-TP-2008- 1160 


Sigtem  Technology,  Inc. 

1343  Parrott  Drive 

San  Mateo,  CA  94402-3630 


6.  AUTHOR(S) 

Chun  Yang  (Sigtem  Technology,  Inc.) 

Mikel  Miller  and  Thao  Nguyen  (AFRL/RYRN) 
Erik  Blasch  (AFRL/RYAA) 


3.  DATES  COVERED  (From  -  To) 

08  April  2005  -  08  April  2007 


1.  REPORT  DATE  (DD-MM-YY) 

April  2007 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

OF  ABSTRACT: 

SAR 

OF  PAGES 

24 

Thao  Nguyen 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

N/A 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


i 


14.  ABSTRACT  (concluded) 


Although  most  of  the  results  are  available  in  the  literature,  they  are  scattered  in  different  textbooks  and  research  papers.  By 
presenting  them  in  a  cohesive  and  concise  manner,  the  paper  can  serve  as  (1)  an  introductory  material  with  key  references  for  those 
who  would  like  to  study  the  subject  further,  (2)  a  quick  reference  for  those  who  would  like  to  compare  these  integration  schemes, 
and  (3)  an  integrated  approach  under  varying  conditions.  Simulation  results  are  presented  to  compare  these  integration  schemes  in 
terms  of  detection  performance  (Pd  and  Pfa)  as  a  function  of  signal-to-noise  ratio  (SNR)  and  integration  interval  under  two 
operating  conditions.  One  is  for  a  signal  with  constant  Doppler  frequency  and  the  other  is  for  a  linear  frequency  modulation  signal 
(i.e.,  a  chirp  signal).  We  show  the  curves  of  the  input  SNR  required  under  different  chirp  rates  so  as  to  maintain  the  same  detection 
performance,  the  resulting  LRTs  for  the  different  approaches,  as  well  as  simulated  comparisons  between  tradeoffs  of  Pd  versus 
SNR. 


11 


Comparative  Study  of  Coherent,  Non-Coherent, 
And  Semi-Coherent  Integration  Schemes 
For  GNSS  Receivers 


Chun  Yang  Mikel  Miller 

Sigtem  Technolog}’,  Inc  Air  Force  Research  Lab/SNRN 

Erik  Blasch  Thao  Nguyen 

Air  Force  Research  Lab/SNAA  Air  Force  Research  Lab/SNRN 


BIOGRAPHIES 

Dr.  Chun  Yang  received  his  title  of  Docteur  en  Science 
from  the  Universite  de  Paris  (No.  XI,  Orsay),  France,  in 
1989.  After  two  years  of  post-doctoral  research  at  the 
University  of  Connecticut,  he  moved  on  with  his 
industrial  R&D  career.  Since  1993,  he  has  been  with 
Sigtem  Technology,  Inc.  and  has  been  working  on 
numerous  GPS,  integrated  inertial,  and  adaptive  array 
related  projects.  He  is  the  co-inventor  of  seven  issued  or 
pending  U.S.  patents.  He  is  also  an  Adjunct  Professor  of 
Electrical  and  Computer  Engineering  at  Miami  University 
of  Ohio. 

Dr.  Mikel  Miller  is  the  Technical  Advisor  for  the 
Reference  Systems  Branch,  Sensors  Directorate,  Air 
Force  Research  Laboratory,  WPAFB,  OH.  He  received 
his  Ph.D.  in  Electrical  Engineering  from  the  Air  Force 
Institute  of  Technology  (AFIT),  WPAFB,  Ohio,  in  1998. 
Since  1986,  he  has  focused  on  navigation  system  R&D 
related  to  GPS,  GPS/INS  integrations,  alternative 
navigation  techniques  including  bio-inspired  navigation 
and  signals  of  opportunity  based  navigation,  autonomous 
vehicle  navigation  and  control,  and  multi-sensor  fusion. 
He  is  currently  responsible  for  both  in-house  and 
contracted  R&D  projects  advancing  navigation 
technology.  He  is  also  an  Adjunct  Professor  of  EE  at 
AFIT  and  Miami  University  of  Ohio.  Dr.  Miller  is  a 
member  of  the  ION,  RIN,  IEEE,  and  AIAA. 

Dr.  Erik  Blasch  is  a  Fusion  Evaluation  Tech  Lead  for  the 
Assessment  &  Integration  Branch,  Sensors  Directorate, 
Air  Force  Research  Laboratory,  WPAFB,  OH.  He  is  an 
adjunct  professor  at  Wright  State  University  (WSU)/ 
University  of  Dayton  (UD)/AFIT  and  a  Reserve  Major  at 
the  Air  Force  Office  of  Scientific  Research  (AFOSR).  He 
received  his  Ph.D.  in  Electrical  Engineering  from  WSU  in 
1999,  a  MSEE  from  WSU  in  1997,  a  MSME  from 
Georgia  Tech  in  1994,  and  a  BSME  from  MIT  in  1992 
among  other  advanced  degrees  in  engineering,  health 


science,  economics,  and  business  administration.  He  is 
the  President  of  the  International  Society  of  Information 
Fusion  (ISIF)  and  active  in  IEEE,  SPIE,  and  ION.  His 
research  interests  include  target  tracking,  sensor  fusion, 
automatic  target  recognition,  biologically-inspired 
robotics,  and  controls. 

Thao  Nguyen  is  a  Ph.D.  candidate  in  the  Electrical  and 
Computer  Engineering  Department  at  the  Air  Force 
Institute  of  Technology  (AFIT),  Wright-Patterson  AFB, 
Ohio.  He  also  works  as  an  Electronics  Engineer  at  the 
Reference  Systems  Branch,  Sensors  Directorate,  Air 
Force  Research  Laboratory,  WPAFB,  OH.  Mr.  Nguyen 
has  been  involved  with  navigation-related  research, 
development,  and  test  since  2002  and  his  current  areas  of 
interest  include  high  anti-jamming  technologies  for  GPS 
receivers,  GPS/INS  integration,  multi-sensor  fusion,  and 
personal  navigation  using  signals  of  opportunity  (SoOP). 

ABSTRACT 

Three  integration  schemes,  namely,  coherent,  non¬ 
coherent,  and  semi-coherent  integration,  and  their  variants 
are  described  in  this  paper  for  signal  detection  and 
tracking  in  such  applications  as  pulsed  Doppler  radar  and 
GPS  receiver.  We  adopt  a  unified  approach  to  present 
such  integration  schemes  as  best  approximations  to  the 
likelihood  ratio  test  (LRT)  when  the  signal  models  are 
partially  known  to  different  extents.  Formulas  for  the 
conditional  probability  density  functions  (pdfs)  of  the 
sufficient  statistic  under  alternative  hypotheses  as  well  as 
the  probability  of  detection  ( Pd )  and  probability  of  false 
alarm  ( Pfa )  are  given.  Integration  improvement  factors  and 
integration  losses  are  then  formulated.  Connections  of  the 
semi-coherent  integration  scheme  to  the  cyclic 
autocorrelation  function,  the  discrete-time  Wigner-Ville 
transform,  and  the  radar  ambiguity  function  are  also 
outlined.  Although  most  of  the  results  are  available  in  the 
literature,  they  are  scattered  in  different  textbooks  and 
research  papers.  By  presenting  them  in  a  cohesive  and 
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concise  manner,  the  paper  can  serve  as  (1)  an  introductory 
material  with  key  references  for  those  who  would  like  to 
study  the  subject  further,  (2)  a  quick  reference  for  those 
who  would  like  to  compare  these  integration  schemes, 
and  (3)  an  integrated  approach  under  varying  conditions. 
Simulation  results  are  presented  to  compare  these 
integration  schemes  in  terms  of  detection  performance  (Pd 
and  Pfa)  as  a  function  of  signal  to  noise  ratio  (SNR)  and 
integration  interval  under  two  operating  conditions.  One 
is  for  a  signal  with  constant  Doppler  frequency  and  the 
other  is  for  a  linear  frequency  modulation  signal  ( i.e .,  a 
chirp  signal).  We  show  the  curves  of  the  input  SNR 
required  under  different  chirp  rates  so  as  to  maintain  the 
same  detection  performance,  the  resulting  LRTs  for  the 
different  approaches,  as  well  as  simulated  comparisons 
between  trade-offs  of  Pd  versus  SNR. 

INTRODUCTION 

In  this  paper,  we  study  several  integration  schemes 
frequently  encountered  for  signal  detection  and  parameter 
estimation  in  such  applications  as  pulsed  Doppler  radar 
and  GPS  receivers.  The  signals  employed  in  these  two 
fields  are  quite  different  in  their  operating  modes  and 
waveforms  (e.g.,  frequencies,  modulations  and  power 
levels).  Even  different  names  are  given  to  similar 
integration  schemes:  pulse  integration  for  radar  vs.  post¬ 
correlation  integration  for  GPS.  Nevertheless,  the 
principles  and  methodologies  are  common  and  our  goal  is 
to  integrate  the  approaches  toward  advanced  detection 
capabilities.  For  an  experienced  reader,  he  or  she  may  go 
directly  to  the  section  Simulation  Results  and  Analysis 
that  highlights  the  comparative  study. 

For  a  signal  buried  in  noise,  a  single  signal  sample  is 
typically  not  suitable  for  detection  simply  because  the 
signal  in  such  raw  samples  is  too  weak  {i.e.,  low  SNR)  at 
a  high  sampling  rate.  To  improve  detection  performance, 
multiple  successive  samples  are  therefore  processed 
jointly,  which  also  serves  the  purpose  of  data 
compression.  In  a  pulsed  Doppler  radar,  a  pulse-matched 
filter  is  applied  to  signal  samples  to  accumulate  the  signal 
energy  per  pulse  while  averaging  out  noise.  Similarly,  in  a 
GPS  receiver,  a  despreading  correlator  is  applied  to  signal 
samples  per  code  epoch,  which  is  1  ms  long  for  GPS  C/A- 
codes.  The  resulting  measurement  per  pulse  or  per  code 
epoch  is  typically  complex-valued,  containing  real  and 
imaginary  parts,  also  known  as  in-phase  and  quadrature- 
phase  components. 

Pulse  measurements  in  radar  or  pre-detection  correlations 
in  GPS  are  further  processed  before  a  decision  is  made  to 
ascertain  the  presence  or  absence  of  a  desired  signal. 
Decision  making  is  carried  out  by  means  of  threshold 
comparison  and  detection  performance  is  characterized  in 
terms  of  probability  of  detection  (Pj)  and  probability  of 
false  alarm  {Pj,,)  for  a  detection  threshold.  One  way  to 
optimally  process  these  measurements  is  to  reduce  them 
into  a  sufficient  statistic,  denoted  by  A.  Under  different 


conditions,  this  leads  to  different  integration  schemes, 
which  are  the  subject  of  study  in  this  paper. 

Indeed,  we  choose  a  unified  way  to  present  various 
measurement-combining  (or  integration)  schemes  as  the 
best  approximation  by  a  sufficient  statistic  of  the 
likelihood  ratio  test  (LRT)  resulting  from  hypothesis 
testing  when  the  underlying  signal  model  is  partially 
known  to  a  different  degree.  There  are  two  hypotheses  of 
interest.  One  is  the  null  hypothesis,  denoted  by  H0 ,  under 
which  measurement  only  contains  noise.  The  other  is  the 
alternative  hypothesis,  denoted  by  Hi,  under  which  the 
signal  of  interest  is  present  in  the  noisy  measurement. 
Once  a  sufficient  statistic  is  given,  a  major  step  that 
follows  is  to  derive  the  conditional  probability  density 
functions  (pdfs)  of  the  sufficient  statistic  under  the  two 
hypotheses,  namely,  P(A\H0)  and  P{A\H{),  respectively, 
from  which  Pd  and  Pfa  are  evaluated  for  a  given  threshold. 
When  closed-form  solutions  are  not  available,  empirical 
formulas  and  numericaPgraphic  means  are  often  used  to 
facilitate  design  and  implementation. 

Coherent  integration  (here  “coherent”  means  matched  in 
time/phase)  works  on  complex-valued  measurements  in 
which  the  signal’s  phases  from  measurement  to 
measurement  are  in  close  alignment,  thus  allowing 
constructive  addition  of  signal  whereas  the  phase  of  noise 
is  random  from  measurement  to  measurement,  resulting  in 
cancellation.  Coherent  integration  can  attain  the  highest 
integration  gain  ideally  (100%  efficiency).  To  achieve  a 
build  up  of  signal  amplitude,  the  phase  relationship 
between  measurements  has  to  be  known  or  estimated, 
which  is,  however,  sometimes  difficult. 

In  non-coherent  integration,  phase  information  is 
discarded.  It  works  on  real-valued  data  samples  by  adding 
up  the  magnitude  (an  envelope  detector)  or  squared 
magnitude  {i.e.,  the  power,  a  squared  law  detector)  or 
sometimes  log-magnitude  of  measurements.  Because  of 
the  nonlinear  operation,  its  analysis  becomes  more 
complicated.  Squaring  effectively  removes  phase 
variations  from  sample  to  sample  (also  unknown  data  bits 
if  it  exists)  but  it  also  squares  the  noise,  allowing  it  to  add 
up.  Noncoherent  integration  is  less  efficient  than  coherent 
integration  but  it  is  often  used  when  there  are  significant 
mismatches  in  signal  phase  due  to  Doppler  frequency 
errors,  unknown  data  modulation,  and  even  receiver  clock 
drifts. 

Both  coherent  and  non-coherent  integrations  schemes  are 
well  documented  in  textbooks  and  widely  used  in 
practical  systems.  A  relatively  new  one  is  the  semi- 
coherent  integration  scheme,  which  also  has  many  other 
names.  Semi-coherent  integration  can  be  defined  as  the 
sum  of  the  products  of  the  input  signal  with  the  conjugate 
of  a  delayed  copy  of  itself.  Since  all  conjugate  products 
have  a  common  phase  argument  (under  constant 
frequency),  they  add  “coherently.”  Each  product 
contributing  the  signal’s  “power”  to  the  sum,  it  is 
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therefore  termed  “semi-coherent.”  However,  the  delayed 
conjugate  product  operation  also  squares  the  noise,  thus 
allowing  it  to  add.  It  is  less  efficient  than  the  coherent 
integration  but  slightly  more  efficient  than  the  non¬ 
coherent  integration  in  terms  of  integration  loss  and  more 
robust  against  frequency  variations. 

In  practical  implementations,  the  coherent  and  non¬ 
coherent  as  well  as  semi-coherent  integration  schemes  are 
typically  used  in  combination.  That  is,  a  certain  number 
of  coherent  integrations  are  followed  by  another  number 
of  non-coherent  or  semi-coherent  integrations.  The 
division  of  the  total  number  of  measurements  into 
coherent  and  non-coherent  or  semi-coherent  is  a  design 
parameter  that  needs  to  be  trade-studied  for  a  particular 
application  at  hand.  This  aspect  is  not  covered  in  this 
paper.  Another  aspect  that  is  not  covered  either  is  the 
post-detection  integration  (also  called  binary  integration, 
“M  out  of  N”  detection,  or  classifier  fusion).  References 
can  be  found  in  [Richards,  2005;  Kaplan  and  Hegarty, 
2006], 


When  a  single  measurement  is  modeled  as  a  constant 
signal  A  in  a  white  zero-mean  Gaussian  noise  with 
variance  cr,  the  signal  to  noise  ratio  per  measurement  is 
defined  as  ( SNR)\  =  A2/cr  and  it  completely  specifies  the 
conditional  pdfs  and  detection  performance.  When  such  N 
measurements  are  processed  by  coherent  integration,  the 
outcome  is  again  a  constant  signal  NA  in  a  white  zero- 
mean  Gaussian  noise  with  variance  Ncr.  In  this  case,  the 
output  SNR  easily  can  be  calculated  as  ( SNR)a  =  NA 2/ cr 
=  N(SNR)  i .  When  compared  to  the  single  measurement 
SNR  at  the  input,  the  output  SNR  is  higher,  proportional 
to  the  number  of  measurements  integrated,  thus  achieving 
100%  efficiency.  The  ratio  of  the  output  SNR  over  the 
input  SNR  is  called  integration  improvement  factor  in 
radar  [Mahafza  and  Elsherbeni,  2004]  or  processing  gain 
in  GPS,  denoted  by  GCI,  as: 
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For  nonlinear  integration  schemes,  however,  the  output 
SNR  is  not  well  defined.  This  is  because  with  cross- 
product  terms,  it  is  difficult  to  properly  separate  signal 
from  noise  after  such  a  nonlinear  operation  as  squaring. 
As  a  result,  with  noise  amplified  by  signal  in  the  cross- 
product  terms,  nonlinear  integration  schemes  are  less 
effective  and  the  resulting  variance  is  signal-dependent. 

There  are  two  ways  to  quantify  integration  improvement 
for  nonlinear  integration  schemes.  One  intuitive  way  is  to 
define  an  equivalent  SNR  at  the  output  of  a  nonlinear 
integration  scheme  such  that  the  integration  improvement 
factor  can  be  calculated  using  the  same  formula  as  in  Eq. 
(1).  A  rationale  behind  this  approach  is  this.  According  to 


the  central  limit  theorem,  the  sum  of  a  large  number  of 
independent  identical  distributed  (i.i.d.)  samples 
converges  to  a  Gaussian  distribution  with  suitable  mean 
and  variance,  which  can  be  used  to  calculate  the 
equivalent  output  SNR. 

A  more  rigorous  way,  however,  is  to  specify  an 
equivalent  SNR  at  the  input  (i.e.,  a  single  measurement 
SNR)  for  a  nonlinear  integration  scheme.  The  SNR 
required  to  achieve  a  specific  Pj  given  a  particular  Pfa 
after  N  samples  are  integrated  using  the  nonlinear 
integration  scheme  is  denoted  by  (SNR) NLI,  which  is 
smaller  than  (SNR) i.  They  are  related  by: 

(SNR)nu=(SNR\xGnl[(N)  (2) 

where  Gnli(N)  is  the  nonlinear  integration  improvement 
factor. 

The  nonlinear  integration  improvement  factor  GNLi(N)  for 
a  nonlinear  integration  scheme  is  typically  smaller  than 
the  integration  improvement  factor  GCi(N)  =  N  for 
coherent  integration.  Their  ratio  is  usually  called  the 
integration  loss  defined  as: 

LI  =  GCI  (N)/Gnu(N)  =  N/Gnli(N)  (3) 

Most  of  the  results  of  integration  schemes  are  available  in 
the  literature  but  scattered  in  different  textbooks  and 
research  papers.  In  this  paper,  we  compile  them  within  a 
unified  framework  and  present  it  in  a  cohesive  and 
concise  manner.  We  hope  the  paper  can  serve  as  an 
introductory  material  with  key  references  listed  for  those 
who  would  like  to  study  the  subject  further  and  a  quick 
reference  for  those  who  would  like  to  use  it  in  their  work. 

The  paper  is  organized  as  follows.  We  start  with  a  review 
of  the  major  results  for  detection  of  signal  in  noise  with  a 
single  measurement.  We  then  present  coherent,  non¬ 
coherent,  and  semi-coherent  integration  schemes  as  ways 
to  combining  multiple  measurements,  respectively.  We 
next  analyze  results  of  comparative  simulation  study  and 
finally  conclude  with  a  summary  and  future  work. 

REVIEW  OF  DETECTION  THEORY  WITH 
SINGLE  MEASUREMENTS 

Consider  a  complex  measurement  z  consisting  of  a 
complex  signal  s'  and  a  complex  noise  w.  We  want  to 
determine  if  the  measurement  z  contains  only  noise  w 
(denoted  as  the  null  hypothesis  H0)  or  both  the  desired 
signal  s  and  noise  w  (denoted  as  the  alternative  hypothesis 

HO: 

H0:  z  =  w  (4a) 

H, :  z  =  s  +  w  (4b) 

The  signal  is  modeled  as  s  =  Ad6  with  A  being  a  constant 
amplitude  and  6  a  constant  phase  and  the  noise  w  is  a 
complex  white  Gaussian  noise  with  its  real  and  imaginary 
components  being  of  zero  mean  and  variance  cr. 
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The  probability  density  functions  (pdfs)  that  describe  the 
measurement  to  be  tested  under  each  of  the  two 
hypotheses  are  denoted  by  p(z\H0)  and  p{z\Hx), 
respectively.  For  an  “optimal”  choice  between  our  two 
hypotheses,  there  are  several  criteria  to  choose.  This 
includes  the  Neyman-Pearson  criterion  in  which  the 
decision  is  made  to  maximize  the  probability  of  detection 
Pd  under  the  constraint  that  the  probability  of  false  alarm 
Pfa  does  not  exceeds  a  predefined  value,  the  Bayes 
minimum  cost  criterion,  and  maximization  of  the 
probability  of  a  correct  decision  among  others  [Therrien, 
1989;  Kay,  1998].  The  optimization  under  any  of  these 
criteria  leads  to  the  ubiquitous  likelihood  ratio  test  (LRT): 


&T  A  V^c 

N  N 

X  X 

II 

< 

(5a) 

H\ 

lnA(z)  =  In  piz\H'  ) 
P(z\H0) 

>  In  77 

< 

H0 

(5b) 

Hi 

Mz)  >  Td 

(5c) 
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where  77,  In  rj  and  Td  are  detection  thresholds  for  the 
corresponding  formulations. 

In  Eq.  (5a),  A(z)  is  the  likelihood  ratio  test  and  in  Eq. 
(5b),  InA(z)  is  called  the  log  likelihood  ratio  test.  The 
latter  is  computationally  more  efficient  as  most  pdfs  are  in 
the  form  of  exponentials  and  as  a  monotonic  increasing 
operation,  the  taking  of  logarithm  does  not  affect  the 
detection  performance  (P,j  and  Pfa). 

In  Eq.  (5c),  A(z)  is  called  the  sufficient  statistic,  which,  if 
it  exists,  is  a  function  of  the  measurement  z  and  has  the 
property  that  the  measurement  z  appears  in  the  likelihood 
ratio  or  the  log  likelihood  ratio  only  through  Mz).  In  other 
words,  knowing  the  sufficient  statistic  A(z)  is  as  good  as 
knowing  the  actual  data  z  in  making  an  optimal  decision. 
As  shown  later  in  this  paper,  various  integration  schemes 
are  such  a  sufficient  statistic  under  different  signal 
models. 

The  optimal  detection  problem  now  becomes  solving  for 
the  detection  threshold  from  equations  of  Pd  and  Pt-a 
expressed  in  terms  of  the  sufficient  statistic  X(z )  and 
signal  and  noise  model  parameters.  This  typically 
involves  two  steps.  A  first  step  is  to  manipulate  the 
conditional  measurement  pdfs,  p(z\H0)  and  p(z\Hx),  so  as 
to  obtain  the  sufficient  statistic  pdfs,  denoted  by  p(A.\H0) 
and  p(A\Hx),  respectively.  A  second  step  is  to  evaluate  Pd 
and  Pfa  of  the  sufficient  statistic  as  a  function  of  the 
threshold  Td. 

+co 

Pl  =  \p(X\H])dX  (6a) 

Td 


+00 

pja  =  \Pa\H())ca 

Td 


(6b) 


Ideally,  the  threshold  Td  can  be  “solved”  from  Eq.  (6b)  for 
a  given  level  of  false  alarm  P/a.  The  probability  of 
detection  Pd  is  then  calculated  from  Eq.  (6a)  as  a  function 
of  the  probability  of  false  alarm  Pfa  with  the  threshold  Td 
as  an  intermediate.  This  Pd  vs.  Pfa  relationship  is  the 
receiver  operating  characteristic  (ROC)  curve.  Although 
the  integrals  can  be  expressed  in  terms  of  pdfs,  there  are 
typically  no  closed-form  solutions  except  for  simple 
cases.  Numerical  methods  are  popular  means  used  to 
present  the  results  in  tables  or  plots. 

From  Eq.  (4),  the  measurement  pdf  can  be  written  as: 


p(z)- 


TIG 


G 


(7) 


where  the  superscript  *  stands  for  complex  conjugate. 
Under  H0,  .?  =  0  and  under  Hu  ^  0.  The  log  likelihood 
ratio  is: 


lnA(z)  =  ^y(2Re{s*z}-s*5)  (8) 

cr~ 

where  Re{z}  means  the  real  part  of  z.  Since  the  second 
term  on  the  right  hand  side  of  Eq.  (8)  does  not  depend  on 
the  measurement,  a  sufficient  statistic  can  be  taken  as: 

k(z)  =  Rej/z}  (9) 


Under  Hi  with  .?  =  Aexp(j&),  Eq.  (9)  can  be  further  written 
as: 


k(z)  =  Re]/.?  +  .s*w}  =  A2  +  Re{wAe  j0} 

=  E  +  Re{wAe-je}  (10) 

where  E  =  A2  is  the  total  energy  in  the  expected  signal  s. 

From  Eq.  (9),  let  £=  s*z,  which  is  a  complex  Gaussian.  It 
is  easy  to  show  that  under  H0,  <^~  PM{0,  2  crE}.  Similarly, 
under  Hh  ^  ~  PM  {E,  2 crE}.  Since  the  power  of  the 
complex  Gaussian  noise  splits  evenly  between  the  real 
and  imaginary  parts  of  C  it  follows  that  the  sufficient 
statistic  k  ~  PN~{0,  crE}  under  H0  and  k  ~  PM  {E,  crE} 
under  H }.  In  other  words,  only  half  of  the  noise  power 
competes  with  the  signal  energy  in  this  case.  The 
following  equations  can  be  established: 


P,,=- 


1 _  erf  ( 


-JcTe’ 


pd  =-erfc{erfc-\2Pfa)- 


=  X-erfc{erjc-'(2PJa)  -J(SNR\} 

Td  =erf-'(\-2  Pf,)jME 


(11a) 


(lib) 
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where  (SNR) i  =  A2  hr  =  E/o2  is  the  signal  to  noise  ratio 
for  a  single  measurement,  erf(x)  is  the  error  function,  and 
erfc(x)  =  1  -  erf(x)  is  the  complementary  error  function 
defined  as: 


erf  (x)  =  -2=  [  e  rdt  (12a) 

erfc(x)  =  -f=  f  e~‘  dt  =  1  -  erf  (x)  (12b) 

-W7T  J 


In  the  above  derivations,  it  is  implied  that  the  signal  s  is 
perfectly  known.  However,  it  is  merely  the  case  in 
practice.  It  is  rather  more  reasonable  to  assume  that  the 
initial  phase  is  unknown  and  to  model  the  signal  as 
s  =  3V9  where  the  phase  angle  6  is  a  random  variable 
uniformly  distributed  over  [0,  2  k)  and  independent  of  the 
signal  components  { S' } .  The  pdf  of  H0  is  unchanged  but 
that  of  H\  now  explicitly  depends  on  6 \  denoted  by  p(z\ 
Hu  O)-. 


r  \  V  a\  1  f  (z-seJ  )  (z-seJ  ) 

p(z  I  Hx,0)  =  — vexP( - 2 - -} 

kct  cr 


1  .  z*z  -2 1  f'z  I  cost?  +  's"s ,  nn 

=  — ^exp{ - - '2 - }  (i:S? 
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Using  the  Bayesian  rule  for  random  parameters,  the 
unconditional  pdf  p(z\  Hi)  can  be  obtained  by  averaging 
the  conditional  pdf p(z\  Hh  0)  over  9  as: 


1  2 * 

p(z\Hl)  =  —  \p(z\Hl,0)dO 


1  ZZ  +  JJ,  ,2|j*z| 

=  - Iexp{ - 2 - Ko(  2  ) 


(14) 


where  the  last  term  If')  is  the  modified  Bessel  function  of 
the  first  kind  (zero-order)  defined  as: 

i  In 

I0(z)  =  — —  J  ezays(e)dO 
2k  „ 


The  log-LRT  can  be  written  as: 


In  A  =  In 


2  I  s  z  I 


■\ 
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(15) 


To  avoid  direct  calculation  of  the  Bessel  function  as  well 
as  the  natural  logarithm,  it  is  desired  to  approximate 
ln[/0(")]  with  a  simple  monotonically  increasing  function. 
It  was  shown  in  [Richards,  2005]  that  for  small  x,  ln[/0(x)] 
~  x2/4  and  for  large  x,  ln[/0(x)]  ~  x.  In  other  words,  the 
optimal  detector  is  well  approximated  by  a  magnitude 
squaring  operation  or  a  square  law  detector  for  x  <  5  dB 
whereas  by  an  envelope  detector  when  x  >  10  dB.  A 
practical  sufficient  statistic  can  thus  be  defined  as: 


Envelope  Detector: 

H, 


X(z )  =  |  j  *z  | 


> 

< 

H„ 


r  ,  |  ?*z  |  >  10  dB 

±  d 


(16a) 


Square  Law  Detector: 

X(z)  =  |  s' z  \2  >  Td  ’  I  ?*z  |  <  5  dB  (16b) 

< 

Ho 

The  performance  of  the  envelope  detector  specified  in  Eq. 
(16a)  can  be  established  by  noting  that  ?*z  ~  Jf{0, 
2 crE}  under  Hn  and  ~  3f  {A2,  2 crE)  under  Ht.  Then, 
under  H„,  X  =|  s'z  \  is  Rayleigh  distributed: 


p(X\H0)- 
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And  under///,  X  is  Rician  distributed: 


P(X\H1)  =  \e<j1 
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From  Eq.  (17),  the  probability  of  false  alarm  can  be 
calculated  as: 

+00 

PjJX  \H0)=  \p(X\H0)dX  =  exp 

Td 

Inverting  this  equation  gives  the  detection  threshold  as: 

Td  =  ^~  Ea1 1  n  Pfa  (20) 


Z7! 

Eg1 


(19) 


From  Eq.  (18),  the  probability  of  detection  can  be 
calculated,  with  the  threshold  from  Eq.  (18),  as: 


PAX  \ho  =  Q,\. 


v 


=  Q,lQ(SM)l,J-\nPfa) 


(21a) 


where  Qui',  ■)  is  the  Marcum’s  Q  function  defined  as: 


+oo 

Qm(ci,t)  =  Jtexp[-|(t2  +a1y}a(at)dt 


(21b) 


Compared  to  the  coherent  detector  given  in  Eq.  (9),  the 
envelope  detector  requires  higher  SNR  to  achieve  the 
same  performance.  This  extra  SNR  required  to  maintain 
the  same  detection  performance  is  called  the  detector  loss, 
a  price  paid  for  not  knowing  the  signal  phase.  With  Pfa  = 
10"6,  the  envelope  detector  requires  0.6  dB  higher  SNR 
than  the  coherent  detector  to  achieve  Pj  =  0.9  and  about 
0.7  dB  more  atPd=  0.5. 

Instead  of  averaging  the  Hfs  pdf  over  unknown 
parameters  to  obtain  an  unconditional  pdf  as  in  Eq.  (14), 
another  approach  is  to  use  the  generalized  likelihood  ratio 
test  (GLRT),  in  which  the  likelihood  ratio  is  expressed  as 
a  function  of  the  unknown  signal  parameter.  In  addition  to 
averaging  the  GLRT  over  the  region  of  the  unknown 
parameter,  which  may  be  assumed  uniformly  distributed 
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if  there  is  no  better  model  available,  another  approach  is 
to  obtain  an  estimate  of  the  parameter  by  maximizing  the 
likelihood  ratio  over  the  unknown  parameter  and  then  to 
compare  the  maximized  likelihood  ratio  to  the  threshold 
for  detection.  The  maximization  can  be  done  analytically 
or  numerically  and  the  latter  is  the  theoretic  basis  for 
parametric  search  (sequential  or  parallel)  in  detection. 

Another  signal  parameter  that  is  usually  unknown  is  the 
signal  amplitude  A,  which  is  assumed  to  be  a 
deterministic  scale  factor.  Eqs.  (lib)  and  (21a)  do  not 
depend  explicitly  on  A  (or  the  noise  power  cr)  but  rather 
on  their  ratio  (SNR) i  for  evaluating  the  detection 
performance.  However,  a  practical  implementation 
requires  specific  values  to  set  up  the  detection  threshold. 
One  is  to  use  a  reference  signal  with  unit  energy  s  such 
that  .9  =  seie  =  AselS  .  For  weak  signals,  another  approach 
is  to  implement  a  locally  optimum  detector  (LOD)  for 
unknown  A,  which  maximizes  the  slope  of  the  likelihood 
ratio  at  A  —  0  while  keeping  a  fixed  Pta  [Kassam,  1988]. 


COHERENT  INTEGRATION  -  OPTIMAL 
COMBINING  OF  MULTIPLE  MEASUREMENTS 

After  introducing  the  signal  model,  this  section  first 
formulates  the  perfect  case  with  known  signal.  It  then 
presents  the  formulation  with  unknown  signal  phase. 
Direct  summation  is  then  presented  followed  by  the  use  of 
FFT  as  an  approach  to  coherent  integration. 

Signal  Model 

Consider  N  measurements  where  each  measurement  zn, 
made  at  time  index  /?,  has  the  following  model: 

zn=sn  +  wn ’  n  =  Q,  (22) 


where  sn  is  the  signal  and  vv„  is  a  complex  white  Gaussian 
noise  with  Re{w„}  ~  ON(0,  cr)  and  Im{w„}  ~  ON(0,  cr). 
Putting  the  N  measurements  into  a  vector  format  gives: 

z  =  [z0,  zN_J  (23a) 

*  =  fo»  •••>  sn-J  (23b) 

w  =  [w0,  wNJ  (23c) 

The  joint  pdf  for  N  complex  measurement  samples  is: 


p(z)  =  - 


1 


-exp{- 


■(z-s)H(z-s). 


(24) 


where  H  is  the  Hermitian  (conjugate  transpose)  operator. 
General  Formulation  for  Perfect  Signal  Model 
The  log  likelihood  ratio  is: 

(25) 


lnA(z)  =  —  (2Re{,s;ffz}  -sH  s) 
o~ 


A  sufficient  statistic  can  be  defined  as: 


A(z)  =  Re]^  z} 


(26) 


From  Eq.  (26),  it  can  be  seen  that  the  LRT  leads  to  a 
rather  simple  sufficient  statistic,  which  specifies  the  data 
processing  to  be  performed  on  the  measurements.  It 
affords  several  interpretations  [Richards,  2005].  The  term 
,v"z  is  a  dot  product  of  the  complex  vectors  s  and  z.  It  can 
be  viewed  as  a  correlation  operation.  This  dot  product 
represents  an  FIR  filtering  operation,  evaluated  at  the 
particular  instant  when  the  equivalent  finite  impulse 
response  s  completely  overlaps  with  the  data  vector  z.  It 
can  also  be  viewed  as  a  matched  filter,  thus  detecting  the 
presence  of  the  mean  vector  in  the  data. 

Under  H\  with  sn  =  A, ,exp(/' </>„),  Eq.  (26)  can  be  further 
written  as: 

A(z)  =  Re  {sH  s  +  sH  w}  =  Re{§4.l2  +  ZWM} 

n= 0  n= 0 

=  E  +  Re^wnAne-]*"}  (27) 

n= 0 

N- 1 

where  |2  is  the  total  energy  in  the  expected 

n= 0 

signal  s  and  E  =  NA 2  in  the  equal  mean  case  with  An  =  A. 

Let  £"  =  s;Hz,  which  is  a  complex  Gaussian  because  it  is  a 
sum  of  Gaussian  random  variables.  It  is  easy  to  show  that 
under  H0,  £  ~  27Sf{0,  2 crE}.  Similarly,  under  Hh  £  ~  ON 
{E,  2 crE}.  The  sufficient  statistic  A  ~  ON{0,  crE }  under 
H0  and  A  ~  ON {E,  crE }  under  H\.  Eqs.  (11a),  (lib),  and 
(11c)  for  the  probability  of  false  alarm,  probability  of 
detection,  and  detection  threshold  respectively,  are 
directly  applicable. 

Eq.  (26)  is  called  a  coherent  detector  since  it  assumes  the 
perfect  knowledge  of  the  signal  s.  From  Eq.  (26)  or  Eq. 
(27),  the  signal  power  after  the  operation  (a  dot  product) 
is  A2N 2  whereas  the  noise  variance  is  Ncrl,  ■  The  SNR  at 

the  coherent  detector’s  output,  denoted  by  (SNR)ah  is 
then  related  to  (SNR)  i  by: 

(SNR)CD  =  4-  =  4  =  ^  =  N(SNR\  (28) 

crE  cr"  cr 

This  improves  the  single  measurement  SNR  by  a  factor  of 
N. 

General  Formulation  with  Unknown  Phase 

In  practice,  however,  the  signal  s  is  not  perfectly  known. 
Assume  that  the  initial  phase  is  unknown  and  the  signal 
can  be  written  as  s  =  ‘seiB  where  the  phase  angle  0  is  a 
random  variable  uniformly  distributed  over  [0,  2;r)  and 
independent  of  the  signal  components  { s  } .  The  pdf  of 
H0  is  unchanged  but  that  of  / /  now  explicitly  depends  on 
6,  denoted  by  p(z\  Hu  0): 

p(z  |  HVG)  =  Nl  2N  exp{-  \  (z-seje)"(z-seje)} 

KG  G 
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(29) 


(33) 
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exp{--^-(zffz-2  |  z  |  cos0  +  fHf)} 
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which  can  be  averaged  over  6,  leading  to  the 
unconditional  pdf/?(y|  Hi): 


1  2 * 

p(z\Hl)=  —  \p(z\H],0)d0 


r^expH-rg" z+1HD}io(2  ^  ~2  Z  b  ^30) 
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The  log-LRT  is: 


In  A  =  In 


f2\ sHz'^ 
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(31) 


Similar  to  Eqs.  (16a)  and  (16b),  the  same  approximations 
can  be  used  to  obtain  the  sufficient  statistic  as: 


Envelope  Detector:  _  ~#z  > 

Td 

(32a) 

_  < 

n„ 

Square  Law  Detector:  z  2 

> 

Td 

(32b) 

< 

H0 

The  performance  of  the  envelope  detector  specified  in  Eq. 
(32a)  can  be  established  by  noting  that  z~  /Af{ 0, 
2crE\  under  H0  and  ~  /Af{E,  2 crE)  under  Hi  where  E  is 
the  total  energy  in  J  .  Then,  X  =|  ][H  z  |  follows  a  Rayleigh 
distribution  under  H0  as  given  in  Eq.  (17)  and  a  Rician 
distribution  under//!  as  given  in  Eq.  (18).  Similarly,  Eqs. 
(19),  (20),  and  (21a)  for  the  probability  of  false  alarm, 
detection  threshold,  and  probability  of  detection, 
respectively,  are  directly  applicable  [Skolnik,  2001]. 

Direct  Summation  for  Constant  Signal 

In  the  above  analysis,  the  detection  of  a  signal  in  noise  is 
treated  as  a  likelihood  ratio  test.  The  formulation  of  joint 
detection  with  multiple  measurements  leads  to  the 
structure  of  coherent  integration  (a  dot  product/a  matched 
fdter/a  correlator/an  finite  impulse  response  (FIR)  filter) 
followed  by  detection.  A  coherent  detector  is  used  when 
the  signal  model  is  perfectly  known  whereas  an  envelope 
detector  or  a  square  law  detector  is  used  when  the  initial 
phase  of  the  signal  is  not  known.  In  the  following,  we 
discuss  two  practical  implementations  of  coherent 
integration.  They  are  (1)  direct  summation  and  (2) 
coherent  integration  with  the  fast  Fourier  transform 
(FFT). 

With  the  noise  characteristics  known,  it  is  the  signal 
model  ( i.e .,  its  structure  and  parameters  and  our 
knowledge  of  it)  that  determines  the  solution  of  optimal 
detection.  When  the  signal  is  constant,  the  direct 
summation  of  N  measurements  is  given  by: 


Cds 


2>„=lA= 

71  =  0 


N- 1  N- 1 


^ Sn 

72  =  0  72  =  0 


where  is  an  Nx  1  vector  of  all  ones. 


From  the  above  analysis,  it  is  clear  that  A,(z)  =  Rej^bx}  is 
the  optimal  detector  for  the  case  where  s„=  A.  The  signal 
amplitude  A  does  not  appear  in  the  detector  but  rather  a 
reference  signal  of  unit  is  used  instead.  When  the  phase# 
is  unknown,  the  optimal  detector  can  be  either  X(z)  =  \fDS\ 
or  \{z)  =  \fos\2  depending  on  the  magnitude  of  the  sum. 


For  the  optimal  detector  in  Eq.  (33)  based  on  the  model  of 
a  constant  signal,  we  are  now  to  analyze  its  behavior 
when  the  signal  is  actually  a  sine  wave  with  a  constant 
frequency  fd.  Under  this  condition,  Eq.  (33)  can  be  written 
as: 


Cds  =  AeW{jtfdTfN  -\)} 
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TV-1 

+Xw„ 


72=0 


=  A  exp  {jjfdTs  (A  - 1)} 


sin (rfdTsN)  ,  sin^r;)  y 
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N- 1 


72  =  0 


(34) 


When  the  Doppler  frequency  is  small,  fjTJV  ~  0  and 
sinc(nfdTsN)  =  s in( nfd TsN)/nfd TSN  ~  1.  The  direct  sum  can 
be  simplified  into 

CDS  *  exp {MdTs(N  - 1  )}AN  + 1>„  «  AN  +  f>„  (35) 

72=0  72=0 


The  signal  power  after  the  direct  sum  (only  consider  the 
real  component)  is  A2N2  whereas  the  noise  variance  is 
Na 2  ■  The  direct  sum’s  SNR,  denoted  by  (SNR)ds,  is  then 
given  by  ( SNR)os  =  A2N2/Nc r2  =  N A2  / a2  =  A^SNRh. 
It  improves  the  input  SNR  by  a  factor  of  N.  However, 
when  the  Doppler  frequency  is  large,  the  SNR 
improvement  vanishes  rather  quickly  in  the  form  of  a 
sine-function. 


Coherent  Integration  with  FFT 

For  a  signal  with  a  constant  Doppler  frequency,  the  LRT 
can  be  first  averaged  over  the  unknown  initial  phase  and 
then  expressed  as  a  function  of  the  unknown  frequency  fd. 
Averaging  again  over  fd  within  the  unambiguous  interval 
from  0  to  1  /Ts  (or  -1/27],  to  1/27],)  leads  to  the  following 
sufficient  statistic: 


l/j1, 

Ca(s)=  J/„(2  | -zH l(f)\)df 
0 


(36a) 


where  J  is  a  signal  replica  defined  as: 

Z(/)  =  [So(A  %-,(./ 'if  (36b) 

(L )  =  exp  {j2jfkTsn }  (36c) 

with  //,  being  a  test  frequency  (which  also  differs  from  the 
true  signal  s  in  amplitude  and  initial  phase). 
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For  practical  implementation,  the  integral  in  Eq.  (36a)  is 
replaced  with  a  filter  bank  with  N  filters  as: 


p(r,<p)  =  - - ^exp{- 

2  7ZO 


2a1 


:>exp 

(7 


(40) 


CcM)  =  ^'YJh(2\zHi(fk)\), 

N  k=0 

with  f  =  Jj_,k  =  0, N-\  (37) 

NTs 

The  test  function  can  be  approximated  by  the  following 
sufficient  statistic: 

Aa(z)  *  max  |  zCI(ft)  |*  max  |  s_{fk)H  z  |  (38) 

k  k 

The  coherent  sum  in  Eq.  (38)  can  be  viewed  as  successive 
rotation  of  the  phase  of  measurements  by  an  amount 
specified  by  the  signal  replica  so  as  to  align  the 
measurements  to  the  same  phase  prior  to  summation.  It  in 
fact  implements  a  matched  filter  with  a  filter  parameter  fk. 
Since  the  frequency  is  unknown,  it  is  searched  over  a 
range  of  possible  values  either  in  parallel  or  in  sequence. 

When  the  filter  parameter  is  discretized  into  a  set  of  N 
parameters  with  a  frequency  spacing  inversely 
proportional  to  the  signal  interval,  i.e.,fk  =  k/NTs,  k  =  0, 
...,  N- 1,  it  represents  a  conventional  filter  bank,  which 
can  be  realized  with  FFT.  For  the  filter  in  which  fk  ~  fa,  its 
output  is  close  to  Eq.  (32)  with  a  SNR  increased  by  a 
factor  of  proportional  to  N. 

NON-COHERENT  INTEGRATION  SCHEMES 

The  amplitude  and  phase  of  a  signal  to  be  detected  are 
typically  unknown.  Ignoring  the  amplitude  in  the  LRT 
does  not  affect  the  detection  performance  and  a  reference 
signal  of  unit  amplitude  is  typically  used.  The  unknown 
phase  can  be  averaged  out  from  the  conditional  pdf  or 
from  the  LRT,  leading  to  an  envelope  detector  or  a  square 
law  detector  as  in  Eq.  (32).  Intuitively,  taking  magnitude 
(or  squared  magnitude  or  log-magnitude)  effectively 
removes  the  unknown  phase  (as  well  as  phase  modulation 
if  it  exists  on  the  signal).  After  the  phase  information  is 
discarded,  combining  magnitude  samples  is  referred  as 
non-coherent  integration,  which  is  analyzed  below. 

The  measurement  model  in  Eq.  (22)  can  be  expressed  in 
terms  of  its  real  and  imaginary  parts  or  equivalently  in 


magnitude  and  phase  as: 

z„  =  s„  +  w„  =  4  +  jl„  =  r«  exp  {j<P„}  (39a) 

where 

4  =  Re{^„  +  w„)  =  r  cos(p„)  (39b) 

1,,  = Im isn  +w„}  =  r„  sin(t z>„ )  (39c) 

rn  =  \zn\  =  y](Re{su  +  w„ })2  +  (Im{s„  +  wn})2  (39d) 

<p„  =  tarf 1  (Im  +  w„ }  /  Re  {s'n  +  wn } )  (39e) 

r2  =  z  z  =eW  (39f) 

n  n  n  z>n  In  v  y 

The  joint  pdf  for  r  and  (p  is  given  by: 


The  pdf  for  r  alone  is  obtained  by  integrating  Eq.  (40) 
over  (p,  yielding: 

P(>j  =  \p(r,<p)d<p  =  —rexp{-  ’  -}/0(— t)  (41) 

•  a  la  a~ 


Similarly,  the  pdf  for  cp  alone  can  be  obtained  by 
integrating  Eq.  (40)  over  r. 

Under  H0,  A  =  0  (noise  alone),  Eq.  (40)  becomes  the 
Rayleigh  probability  density  function: 


p(r)=  \  exp{-  ’’  J 
( 7  A- (7 

(42a) 

p((p)  =  —  >  (p  e  (0,  2k) 

2  K 

(42b) 

Under//!,  A^O,  Eq.  (41)  is  the  Rician  probability  density 
function.  Also  when  A/cr  is  very  large,  Eq.  (41)  becomes 
a  Gaussian  probability  density  function  with  mean  A  and 
variance  cr. 


Let  r  =  [/'o,  ...,  rw_!]rbe  a  vector  of  A  magnitude  samples. 
The  joint  pdfs  under  H0  and  H\  are 


P(l  I  R»)  =  n-TexP'_TT} 

«=0  ^ 

y-1  y  y2  +  A2  y  A 

P(r\  A,)  =  n ~^-exp{-  "  2  }/0(~V) 

n= 0  &  2.C 7  (7 


(43a) 

(43b) 


With  the  same  approximation  as  used  to  derive  Eq.  (32b) 
or  (16b),  the  log-LRT  becomes 


InA  =  ^ln 


r/-4) 

„=0  L  O’ 


ir  A 


a 


(44) 


Combining  all  constants  into  the  threshold  leads  to  a 
practical  sufficient  statistic  as: 


^=Yj 


N- 1 

Zl  A  i2 


n= 0 


< 
H  o 


(45) 


which  specifies  a  noncoherent  integration  detection  rule 
using  the  square  law  detector.  A  similar  noncoherent 
integration  detection  rule  using  the  envelope  detector  can 
result  if  the  log-LRT  is  approximated  using  Eq.  (32a)  or 
(16a). 


To  determine  the  detection  performance,  we  start  with 
calculating  the  pdfs  for  the  sufficient  statistic  in  Eq.  (45). 
First,  let  q  =  r2  and  we  determine  the  pdfs  for  q  from  Eq. 
(43)  as: 


p(q) 


l 

2c? 


exp{- 


q  +  A 2 
2cr 


-KoO 


V?4 

2  ' 


(46) 
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When  A  =  0,  q  =  r2  =^2  +  tj2  obeys  the  Chi-square  (yj) 
distribution  with  two  degrees  of  freedom,  which  is  also  an 
exponential  distribution.  For  A  ^  0,  it  obeys  a  non-central 
Chi-square  (%  )  distribution.  When  the  degrees  of  freedom 
exceed  30,  a  non-central  Chi-square  (y2)  distribution  can 
be  accurately  approximated  by  a  Gaussian  distribution. 

Since  the  random  variables  r„  (and  the  corresponding  qn  = 
r  )  for  all  /?  are  independent,  the  pdf  for  the  sufficient 

statistic  A  can  be  derived  as  the  convolution  of  individual 
pdfs: 

A  —  A 

!  exp,-s7-JV(S™)^ 

<«) 


Non-coherent  integration  is  less  efficient  than  coherent 
integration.  The  non-coherent  integration  gain  is  always 
smaller  than  that  of  the  coherent  integration  for  the  same 
number  of  measurements.  This  integration  loss,  estimated 
somewhere  J~N  between  and  A,  is  approximated  by  a 
simple  formula  [DiFranco  and  Rubin,  1980]: 

LNCI  =  101og(Viv")-5.5 ,  dB  (52) 


Another  approximation  for  the  integration  loss  factor 
[Barton,  1988;  Curry,  2001]  is: 


4vCT  ~ 


1  +  (SNR), 
C SNR), 


(53) 


Yet  another  approximation,  which  is  accurate  within  0.8 
dB  given  by  [Peebles,  1998],  is  in  terms  of  integration 
improvement  factor  as: 


where  ( SNR),  =  A2!cf  as  previously  defined  and  /,v-i C-)  is 
the  modified  Bessel  function  of  order  A- 1. 

When  A  =  0,  the  pdf  for  the  sufficient  statistic  A  is  yj- 
distributed  with  2 A  degrees  of  freedom,  which  can  be 
written  as: 


p(A\H0)  = 


1 


c t2N2nT{N ) 


AN~'  exp{-— 
1(7 


(48) 


where  T\N)  =  (A- 1 ) !  is  the  Gamma  function.  As  a  special 
case  of  the  gamma  density,  the  Erlang  density  reduces  to 
the  exponential  pdf  when  A  =  1. 

The  probability  of  false  alarm  is  obtained  by  integrating 
the  pdf  in  Eq.  (48)  from  the  threshold  to  infinity  as: 


Pfa(A \H„)  =  J p(A  |  H„)dA  =  1 

Td 


(49) 


where  the  last  term  is  called  Pearson’s  form  of  the 
incomplete  gamma  function  defined  as 


l{u,M) 


J 


Ml 


dr 


(50) 


G(N)nci  =  6.79(1 +  0.235PJ 


1  + 


logjl/i^)' 

46.6 


logic  (A)  x 


[1  -  0.141og10(A)  +  0.018311og[0(A)] ,  dB  (54) 


Using  Eq.  (53),  the  output  SNR  after  non-coherent 
integration  can  be  written  as: 


(SNR)nci  = 


( SNR)a 

4vO 


=  N(SNR\ 


_  N(SNR)l 

L NCI 

(SNR), 

1  +  (SNR), 


(55) 


Strictly  speaking,  the  integration  loss  factor  or 
improvement  factor  ought  to  be  specified  by  calculating 
an  equivalent  SNR  at  the  input  (i.e.,  a  single  measurement 
SNR)  that  is  required  to  achieve  a  specific  Pj  given  a 
particular  Pfa  after  A  samples  are  integrated  using  a 
nonlinear  integration  scheme.  This  is  because  the  output 
SNR  for  a  nonlinear  integration  scheme  is  not  well 
defined  due  to  cross-product  terms. 

Flowever,  some  equivalent  SNR  is  frequently  used  at  the 
output  of  a  nonlinear  integration  scheme.  One  definition 
is  given  in  [Kay,  1998]  as: 


The  probability  of  detection  is  obtained  by  integrating  the 
pdf  in  Eq.  (47)  from  the  threshold  to  infinity  as: 


Pd(A\H^=)p(A\H0)dA  =  Qk 


2NA 


<7 


-(r+^L)  n 


+  e 


I 


NA2 


A„-i  (2 


NTdA* 

a 2 


)  (51) 


Closed-form  solutions  generally  do  not  exist  for  these 
integrals.  Numerical  solutions  are  often  used  to  generate 
tables  and  curves.  Flowever,  some  empirical 
approximations  exist  for  manual  calculation.  One 
example  is  Albersheim’s  equations  [Richards,  2005]. 


E{z\H,}-E{z\H0}  (56) 

JVar{z\H0} 


where  the  square  root  is  used  for  variance  in  the 

denominator  so  that  the  output  SNR  is  the  ratio  of  the 

same  units  as  the  input  SNR. 

As  another  example,  consider  the  non-coherent 
integration  scheme  using  an  energy  detector  in  Eq.  (45).  It 
is  easy  to  show  that 

E{A\H0}  =  2Na2  (57a) 

Var{A  \H0}  =  8Acr4  (57b) 

E{A\H,}=N(A2  +2cj2)  (57c) 
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According  to  Eq.  (56),  the  SNR  for  the  sufficient  statistic 
in  Eq.  (45)  is: 


(SNR)nci_fd  = 


A2 

2cr 


(58) 


The  following  approximations  are  given  by  [Lowe,  1999]: 


(SNR)n 


4-x(  4 


64 


768 


(64a) 


Compared  to  coherent  integration,  the  corresponding 
integration  loss  is: 

Lnci  ED=JlN  «\a4n  (59) 


which  is  consistent  with  the  loss  estimates  given  before. 

Consider  yet  another  example  where  the  integration  is 
done  of  A  outputs  from  an  envelope  detector  as: 


N-\ 

El  A  I 


n= 0 


< 

Ho 


(60) 


which  can  be  viewed  as  if  the  log-LRT  Eq.  (44)  is 
approximated  using  Eq.  (32a).  The  envelope  detector’s 
output  follows  a  Rayleigh  distribution  Eq.  (42a)  under  H0 
and  a  Ricean  distribution  Eq.  (41)  under//!,  respectively. 
The  conditional  means  and  variances  can  be  computed  as: 


[x  1  1  3  V 

\2  +  2(SNR),  +  8(SNR)i'z  +  16 {SNR)?2  y 


(iSSVH),  >  2.8  (64b) 


When  (SNR)  i  «  1,  the  signal  is  too  weak  to  detect  from  a 
single  measurement.  The  lower  order  term  in  Eq.  (64a) 
indicates  that  the  output  SNR  is  proportional  to  (SNR)  | 
but  to  the  square  root  of  A  at  this  regime.  From  Eq.  (64a), 
the  loss  can  be  written  as: 


^A  «E05VA 


(65) 


Other  effects  that  further  introduce  SNR  loss  include 
amplitude  fluctuations  such  as  Swelling  models 
[Richards,  2005]  and  phase  noise  [Richards,  2003],  which 
are  not  covered  in  this  paper. 


SEMI-COHERENT  INTEGRATION  SCHEMES 


E{r\Ha}  =  A-a 


Var{r\H0}  =  E{r2\H0}-E2{r\H0}  = 


4 -71 


E{r\Hi}  =  i\-oe 


A2  A2  A2  A2 


Var{r  \HX}  =  A2  +  la2  -  E2  {r  \  Hx } 


(61a) 
(61b) 
(61c) 
(6  Id) 


Both  Eqs.  (61a)  and  (61b)  can  be  derived  from  Eqs.  (61c) 
and  (61d)  by  setting  A  =  0.  Assume  that  \rn\  are  i.i.d.  and 
then  the  conditional  means  and  variances  of  A  can  be 
computed  from  those  of  r„  in  Eqs.  (61a)  through  (61d)  as: 


For  a  signal  z  of  length  A  with  an  amplitude  A,  there  are 
an  unknown  initial  phase  <f b  and  an  unknown  Doppler 
frequency/}.  The  initial  unknown  phase  can  be  averaged 
out  as  done  previously.  The  resulting  likelihood  ratio  can 
be  further  averaged  over  the  phase  increment  caused  by 
the  unknown  Doppler  per  sampling  interval,  denoted  by  (p 
=  2 7tfdTs  (the  total  phase  is  <j>n  =  <fo  +  (p  «),  which  may  be 
assumed  uniformly  distributed  over  [0,  2ji).  The  log-LRT 
can  be  written  as: 


LnA  =  — 
2  n 


\l0(2A 


0 


zH1(<p)\)d(p 


(66) 


E{A\H0}  =  NE{r\H0}  =  NA-cr 


Var{A  |  H0}  =  NVar{r  \  H0}  =  A 


4 -71 


(62a) 

(62b) 


E{A\Hl}=NE{r\Hl}=N^-oe 

(1  +  Eiw  (El)  +  El/  (Ei)l  (62c) 

(  2a:)  o(2a2)  2a2  l(4a2)J 

Var{A  |  Ht)  =  NVar  {r  \  Hx}  =  N(A2  +  la2  -E2{r  \  Hx})  (62d) 

The  output  SNR  of  non-coherent  integration  with  an 
envelope  detector  according  to  Eq.  (60)  is: 


(SNR)nci  ED 

=  4n. 


=  VAJ-^ 

4- n 


A  A  A2  A\ 
(  +  2 A  0  2ct2  +  2cr2  4<r! 


-1 


a+«  vo(« )+« /i(Ml) 

V  2  2  2  4 


-1 


(63) 


The  Bessel  function  I0(m)  in  Eq.  (66)  can  be  approximated 
by  a  Taylor  series  for  small  values  of  A.  According  to 
[Selin,  1965],  the  result  is: 


LnA  = : 


64 


l6  5  . 

A2  0 


-2 


(67) 


where 


N 

sk  =  XZ»ZE 


n=k+ 1 


(68) 


The  only  term  in  Eq.  (67)  that  depends  on  the  signal’s 
Doppler  phase  is  the  sum  with  Sk  for  k  *  0  whereas  the 
other  two  terms  involve  the  total  power  So  from  both  the 
signal  and  noise.  It  was  suggested  in  [Wirth,  2001]  to  use 
the  sum  with  Sk  in  Eq.  (68)  as  a  test  function: 

r,  =  af  « 

k= 1 
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which  is  also  called  the  autocorrelation  estimate  (ACE) 
test  because  Sk  is  related  to  the  ACE  of  the  signal  series 
{ z„j  for  the  time  lag  kTs  by: 


Pk  = 


1 


N-k 


N 

TjZnZl-k 


n=k+ 1 


1 


N-k 


Sk 


(70) 


The  parameter  K  is  at  most  N- 1 .  When  K  =  0,  we  have  k  = 
0  and  T0  =  S0.  This  corresponds  to  non-coherent 
integration  (with  a  square  law  detector  or  by  power- 
combination)  discussed  in  the  previous  section.  K  =  1 
leads  to  k  =  1,  which  corresponds  to  the  so-called  tone 
detector  to  be  discussed  later.  In  general,  the  detection 
performance  improves  with  increasing  K.  However,  with 
a  large  k,  the  number  of  samples  available  for  calculating 
Sk  is  reduced.  Typically,  KIN  =  'A  ~  14. 

For  a  signal  with  constant  Doppler  zn  =  Asyzpijtpn),  the 
product  of  the  input  signal  with  the  conjugate  of  a  delayed 
copy  of  itself,  i.e.,znz\_k  =  A2  exp( j<pk),  has  a  phase 

argument  independent  of  n.  As  a  result,  the  summation  in 
Eq.  (68)  adds  z  z*  ,  “coherently”  for  all  the  measurement 

contributions  over  n.  But  what  each  measurement 
contributes  to  the  sum  is  the  signal  “power.”  It  is  in  this 
sense  that  the  summation  is  called  “semi-coherent.” 
Delayed  squaring,  although  removing  phase  variations 
from  sample  to  sample  (and  also  unknown  data  bits),  also 
squares  the  noise,  thus  allowing  it  to  add  up.  It  is 
therefore  less  effective  than  coherent  integration. 

For  a  small  k,  an  additional  SNR  is  required  to  achieve  the 
same  detection  performance  as  an  ideal  coherent 
integration,  which  is  the  integration  loss  mentioned 
earlier.  As  noted  in  [Wirth,  2001],  the  ACE  test  is 
expected  to  operate  on  a  single  signal  and  a  weak  signal 
may  be  suppressed  by  a  stronger  signal  due  to  the 
nonlinear  nature  of  the  test.  It  was  shown  in  [Wirth,  2001] 
that  for  K  =  N/2,  the  loss  compared  to  a  fdter  bank  is  only 
0.5  dB.  Since  a  filter  bank  may  suffer  from  a  straddling 
loss  about  1  dB,  there  may  even  be  a  gain  of  0.5  dB 
compared  to  a  filter  bank.  Another  advantage  of  the  ACE 
test  is  that  it  is  more  robust  against  frequency  variations. 
In  the  next  section,  we  will  show  these  results  via 
computer  simulations. 

Now  consider  the  special  case  where  K  =  1.  The 
summation  of  z  z*  ,  over  n  for  k  =  1  is: 

n  n-k 


<r=I 


K„Z„ 


(71) 


where  k„  is  a  weighting  factor. 

The  test  function  z  in  Eq.  (67)  is  called  a  tone  detector  in 
[Clarkson,  1997].  The  test  £  for  k„  =  1  is  known  as  the 
semi-coherent  statistic,  which  has  been  studied  since  1960 
[Reed  and  Swelling,  1962;  Lank,  Reed,  and  Pollon, 
1973].  The  statistic  is  also  known  as  the  weighted  linear 


predictor  for  different  choices  of  k„,  [Kay,  1989]  and 
studied  in  [Clarkson,  Kootsookos,  and  Quinn,  1994]. 

Following  the  approach  of  [Clarkson,  1997],  the  additive 
noise  in  Eq.  (39a)  can  be  arranged  into  a  multiplicative 
one  as: 

zn  =exp{/(0  +  2#h)}Q4  +  t>„)  (72) 

where/is  an  unknown  signal  frequency. 

Since  v„  is  related  to  the  original  noise  w„  through 
multiplication  by  an  exponent  with  purely  imaginary 
argument,  it  is  also  a  complex  white  normal  random 
variable  with  the  same  variance. 

The  tone  detector  In  Eq.  (71)  can  be  written  as: 


£  =  e]2,fYjKn(A  +  on+1  )(A  +  u„)  =  eJ2*Z  (73) 

M  =  0 

It  follows  that,  for  Kn  =  1, 

mz=E{Z}  =  A2(N- 1)  (74a) 

<4  =  E\Re(Z -mz)2}  =2a2A2(2N-3)  +  2a\N-l)  (74b) 
(72Zq=E{lm(Z-mz)2}=2a2A2  +  2cr\N-l)  (74c) 

E{Re(Z -mz)Im(Z -mz)}  =0  (74d) 


From  the  above  equations,  it  is  easy  to  see  that  for  A  =  0, 
crZi  =  cr2Zq  =  2 a\N  - 1)  and  a\  =  Aa\N  - 1)  •  According  to 
Eq.  (56),  the  output  SNR  of  the  tone  detector  in  Eq.  (71) 
for  k„  =  1  is: 

A2 

(, SNR)m  =  -2-  4n  =  Jn(SNR),  (75 ) 

2er 

which  is  related  to  the  energy  detector  in  Eq.  (58)  as: 

(SNR)Td  _  k  (76) 

(SNR)nci _ed 

Comparison  of  the  noise-only  variances  shows  that 
a\D  =  4 a4  {N  - 1)  <  <j2ed  =8 a4N  ■  In  other  words, 
° td^ed  =(N-  Y)/2N  ■  This  shows  that  the  SNR  required 
to  achieve  the  same  detection  performance  by  the  energy 
detector  (;'.e.,  the  square  law  detector)  is  yfl  higher  (or 
1.5  dB  higher)  than  that  required  by  the  tone  detector 
[Clarksons,  1997]  and  is  also  the  asymptotic  result  of 
[Lank,  Reed,  and  Pollon,  1973]. 

Recently,  this  kind  of  semi-coherent  integration  scheme 
has  been  used  for  differentially  coherent  PN  code 
acquisition  in  DS-SS  CDMA  receivers  [Zarrabiadeh  and 
Souza,  1997]  and  for  assisted  acquisition  of  weak  GPS 
signals  [Choi  et  al.,  2002].  It  was  estimated  in  [Yu,  2006] 
that  the  improvement  of  sensitivity  with  such  a 
differential  combining  is  1.2  dB  to  1.6  dB  over  the 
noncoherent  integration  for  Pj  =  0.1  and  Pfa  =  10"3,  which 
is  consistent  with  the  above  estimate.  In  fact,  a  more 
rigorous  derivation  based  on  the  joint  probability  density 
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function  of  magnitude  and  phase  for  consecutive 
measurements  was  used  in  [Yu,  2006].  It  also  explained 
by  means  of  Cauchy-Schwartz  inequality  that  the  noise  in 
the  tone  detector  is  always  smaller  than  that  in  the  energy 
detector.  This  can  be  seen  by  comparing  Eq.  (45)  and  Eq. 
(71)  with  A  =  0  as: 

Var{£m  |  H0}  =  E{£m£D  \  H0}  =  E(£  wn+lw„(X  »wO*} 

n= 0  m= 0 

N-2  N- 1 

=  E(£  wn+lwn+lwnwn}  < 

n= 0  n= 0 

V  n= 0  n= 0 

=  £{(£  W>„)2}  =  Var{£EDfED}  (77) 

«=0 

If  the  input  signal  is  a  chirp  signal  (linear  frequency 
modulation),  the  product  terms  in  Eq.  (67)  contains  a 
monotonic  signal  with  a  frequency  of  (2 JSTJc)  where  /?  is 
the  chirping  rate  and  Tsk  is  the  delay  time.  This  is  a  beat 
frequency  between  the  signal  and  its  delayed  copy,  which 
can  be  estimated  by  performing  an  FFT  to  the  sequence  of 
conjugate  products.  The  chirp  rate  can  then  be  determined 
from  the  estimated  beat  frequency  [Zaino  et  al.,  2000]. 
The  test  function  becomes: 


C(m,k)  =  YJznzl-mex Pf 


-  jlTtkn^ 


(78) 


Without  noise,  this  is  equivalent  to  a  matched  chirp  filter. 
The  performance  is  close  to  the  coherent  integration  when 
the  SNR  is  large.  However,  at  low  SNR,  signal 
suppression  can  occur. 

It  is  interesting  to  note  that  the  test  function  specified  in 
Eq.  (78)  is  similar  in  structure  to  the  cyclic 
autocorrelation  function  [Gardner,  Napolitano,  and  Paura, 
2005]  defined  at  frequency  a  for  lag  parameter  x  as: 


Ra{  t)  =  E{x(n  +  ^ )x(n  -  T-)e]lnm} 

=  T7  iLx(n  +  ^)x(n-^)e-J2™  (79) 

JS!  n=\  1  2. 


In  addition,  the  similarity  of  the  cyclic  autocorrelation 
function,  the  discrete  Wigner-Ville  transform,  and  the 
semi-coherent  integration  with  the  radar  ambiguity 
function  may  be  recognized.  However,  the  relationship 
between  them  is  only  superficial.  As  pointed  in  [Gardner, 
1991],  the  concepts,  theory,  and  applications  underlying 
these  definitions  have  little  in  common. 

SIMULATION  RESULTS  AND  ANALYSES 

For  the  simulations  presented  below,  the  following  signal 
model  is  used: 

s„  =  =  A^V{j\2^(fJ,n  +  fiT;n2)  +  $,]}  (81) 

where  A  is  the  signal  amplitude,/)' is  an  unknown  Doppler 
frequency,  /?  is  an  unknown  chirping  rate,  (4)  is  an 
unknown  initial  phase  uniformly  distributed  between  0 
and  2 n,  and  Ts  =  1  lfs  is  the  sampling  interval  with  /  being 
the  sampling  frequency. 

To  simulate  a  desired  SNR  level,  the  signal  amplitude  A  is 
adjusted  according  to: 

C/Wp  _ 

^  =  (7;io  10  y  (82) 

where  C/N0  is  the  corresponding  carrier  to  noise  ratio  in 
dB-Hz. 

In  the  simulations,  a  complex  noise  is  used  in  which  the 
real  and  imaginary  parts  are  real  white  zero-mean 
Gaussian  with  variance  c2  =  l.  Under  the  Nyquist  rate, 
the  signal  passes  through  a  low  pass  filter  with  an 
equivalent  bandwidth  of  ±//2.  For  Ts  =  0.001  s,  A  =  1 
corresponds  to  C/N0  =  30  dB-Hz.  This  is  equivalent  to 
SNR  =  101ogio(H2/cf)  =  0  dB.  Similarly,  A  =  0.5  leads  to 
C/N0  =  24  dB-Hz  and  equivalently  SNR  =  -6  dB.  We  will 
use  the  input  signal  (single  measurement)  SNR  values  to 
plot  simulation  results  whenever  applicable. 

We  generate  the  probability  of  detection  ( Pd )  as  a  function 
of  single  measurement  (input)  SNR  (varying  from  -30  dB 
to  10  dB)  for  a  fixed  probability  of  false  alarm  (P/a  = 
10"3).  Six  integration  and  detection  schemes  are 
compared,  which  are  illustrated  in  Fig.  1.  They  are: 


As  a  final  note,  consider  a  symmetric  window  of  length 
M+ 1  centered  at  time  index  n,  take  a  delay  index  m  that  is 
also  symmetric  relative  to  /?,  and  then  vary  it  over  the  data 
window,  i.e.,  y  y'  for  m  e  [-Ml 2,  M/2].  Apply  the 

Fourier  transform  over  these  delayed  conjugate  products 
over  m,  in  contrast  to  n  in  Eq.  (78),  and  this  leads  to  the 
discrete  Wigner-Ville  distribution  (DWVD)  defined  as: 


WV(n,k)  =  2  X  v„+my;_„,exp{-  j2mn2k ' 


M  + 1 


(80) 


where  the  index  k  corresponds  to  the  frequency  at 
2k!Ts{M+\). 


1  -  Ideal  Coherent:  Aa(z)  =  Re{^zA*} 

«=0 

2  -  Practical  Coherent  with  FFT : 

2c,_fft(z)  =  max{|  FFT {z}  |} 


3  -  Non-Coherent:  = 

n= 0 

N- 1 

4  -  Semi-Coherent  for  First  Lag:  Xsa(z)  =  X1 

n= 0 

5  -  Semi-Coherent  up  to  First  N/2  Lags: 

NI2 

^SCI  _N /  2  (— )  =  Xj 
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Ideal  Coherent 
Integration 


Practical  Coherent 
Integration  with  FFT 


Non-Coherent 

Integration 

Serai-Coherent 
Integration  for 
First  Lag  with  FFT 

Serai-Coherent 
Integration  for 
First  Lag 


Serai-Coherent 
Integration  for 
First  K  Lags 


Fig.  1  -  Block  Diagram  of  Six  Compared  Integration  Schemes 


6  -  Semi-Coherent  for  First  Lag  with  FFT: 

A‘SCI_FFT&  =  max(l  FFT  {dlag  (— l:jV-l  )— 0:tV-2  }  1} 

where  FFT{v }  stands  for  an  FFT  operation  applied  to  the 
vector  v  which  may  be  zero-padded  to  reduce  the  size  of 
frequency  bin  if  necessary,  max{y}  stands  for  finding  the 
maximum  value  of  y,  and  diag(v)  stands  for  a  diagonal 
matrix  with  y  along  the  principal  diagonal. 

We  first  consider  the  case  of  a  constant  Doppler 
frequency  with  /3=  0.  The  signal  model  reduces  to 

sn  =  Acxp{j(27ifdTn  +  </)n)}  (83) 

and  the  product  of  the  signal  and  the  conjugate  of  its 
delay  copy  at  lag  k  is: 

s,/„-k  =  a2  axP{J24dTM  (84) 

which  does  not  depend  on  n  and  the  beat  frequency  is 
therefore  fdTsk  =  -0.25  Hz  for  A  —  1 ,  Ts  =  0.001  s,  and  fd  = 
-250  Hz. 

We  ran  the  six  integration  schemes  with  noise  only 
10,000  times  so  as  to  determine  the  detection  threshold 
that  could  produce  the  desired  probability  of  false  alarm 
at  Pfa  =  10"3  (i.e.,  the  threshold  was  crossed  10  out  of 
10,000  times).  We  then  ran  the  six  integration  schemes 


again  10,000  times  for  each  of  the  SNR  values  with  both 
signal  and  noise  present.  The  probability  of  detection  was 
estimated  as  the  number  of  runs  in  which  the  detection 
threshold  was  crossed,  normalized  by  the  total  number  of 
runs. 

Figs.  2,  3,  and  4  show  the  Pd  vs.  SNR  curves  for  N  =  20, 
40,  and  100,  respectively.  As  expected,  the  curves  are 
shifted  toward  the  left  (toward  low  SNR  values)  by  3  and 
7  dB  for  the  coherent  integration  scheme  as  the 
integration  time  is  increased  by  a  factor  of  2  and  5, 
respectively.  In  these  figures,  the  curves  keep  about  the 
same  shape  and  the  SNR  required  to  achieve  Pd  >  0.5 
maintains  the  following  order: 

Ideal  Coherent  <  Practical  Coherent  with  FFT  <  Semi- 
Coherent  up  to  First  N/2  Lags  <  Semi-Coherent  for  First 
Lag  <  Semi-Coherent  for  First  Lag  with  FFT  <  Non- 
Coherent 

The  practical  coherent  integration  scheme  implemented 
with  FFT  (i.e.,  the  green  curve)  is  behind  the  ideal 
coherent  integration  (i.e.,  the  blue  curve)  by  more  or  less 
2  dB,  that  is,  an  integration  loss  of  2  dB  at  Pd  =  0.9  for 
instance.  The  detection  threshold  of  the  practical  coherent 
with  FFT  is  higher  than  the  ideal  coherent  to  satisfy  the 
same  P/a.  This  is  because  the  optimization  across  all 
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P.  =  0.001,  a  =  2a  =  2,  N  =  20,  runs  =  10000 


Pfa  =  0.001,  a,  =  2a  =  2,  N  =  40,  runs  =  10000 


Fig.  2  -  Pd  vs.  SNR  for  N  =  20  =  0  Hz/s) 


single  measurement  SNR,  10log1Q(A2/a2),  dB 

Fig.  3  -  Pd  vs.  SNR  for  N  =  40  (J3  =  0  Hz/s) 


=  0.001,  a  =  2a  =  2,  N  =  100,  runs  =  10000 


P,  =  0.001,  a  =  2a  =  2,  N  =  100,  runs  =  10000 


1  r 

0.9  - 
0.8  - 
0.7  - 
0.6  - 
0.5  - 
0.4  - 
0.3  - 
0.2  - 
0.1  - 


- coherent 

- non-coherent 

- semi-coherent, 

- semi-coherent, 

k->N/2 
k=1,  fft 

Fig.  4  -Pd  vs.  SNR  for  N=  100  (J3=  0  Hz/s) 


Fig. 


-25  -20  -15  -10  -5  0  5  10 

single  measurement  SNR,  10log10(A2/a2),  dB 

5  —  Pd  vs.  SNR  for  N=  100  (J3=  37.45  Hz/s) 


=  0.001,  a  =  2a  =  2,  N  =  100,  runs  =  10000 


f.  =  200  Hz,  N  =  20,  t  =  1  ms 


Fig.  6 -Pd  vs.  SNR  for  TV  =  100  (/?=  749  Hz/s) 


Fig.  7  -  Loss  of  Integrated  Signal  Peak  for  N~  20 
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frequency  bins  tends  to  maximize  noise.  Besides,  the  ideal 
coherent  only  considers  the  noise  component  in  the  real 
part  whereas  the  practical  coherent  with  FFT  takes  on 
noise  in  both  the  real  and  imaginary  parts  via  an  envelope 
detector.  When  the  signal  frequency  is  between  two 
frequency  bins,  an  additional  saddle  loss  may  occur  even 
though  the  FFT  is  zero-padded  to  four  times  the  original 
size  ( i.e .,  80)  to  improve  the  frequency  resolution  in  the 
simulation.  The  integration  loss  is  a  price  to  pay  for  not 
knowing  the  signal  exactly. 

The  semi-coherent  with  FFT  (the  black  curve)  is  slightly 
better  than  the  non-coherent  (the  red  curve)  but  both  are 
behind  the  semi-coherent  for  first  lag  (the  aqua  curve). 
For  the  case  of  J3  =  0  with  a  constant  Doppler  frequency, 
the  terms  of  conjugate  products  have  a  constant  phase 
term.  As  a  result,  the  semi-coherent  with  FFT  produces 
about  the  same  output  as  the  semi-coherent  for  first  lag 
when  the  signal  is  present.  Flowever,  its  detection 
threshold  is  higher  for  the  same  reasons  explained  above. 
The  resulting  detection  performance  is  about  the  same  as 
the  non-coherent,  which  exhibits  an  integration  loss  of  0.9 
dB  at  Pd  =  0.9  as  compared  to  the  semi-coherent  for  first 
lag. 

The  semi-coherent  up  to  first  N/2  lags  (the  purple  curve) 
is  behind  the  practical  coherent  with  FFT  (the  green 
curve)  with  an  integration  loss  of  about  0.3  dB  for  N=  20 
and  about  0.7  dB  for  N  —  100  and  an  integration  gain  of 
1.2  dB  for  N=  20  and  2.8  dB  for  N=  100  over  the  semi- 
coherent  for  first  lag,  respectively. 

In  the  second  simulation,  we  consider  the  case  with  J3^0. 
This  linear  frequency  modulation  introduces  a  frequency 
variation  and  will  adversely  affect  the  coherent 
integration.  On  the  other  hand,  the  product  of  a  signal  and 
the  conjugate  of  its  delay  copy  can  be  written  as: 

s/n-k  =  A2  exp{j2nfdTsk}  exp {-j0Ts2k2}x 

exp{j2n:2j3Ts2kn}  (85) 

where  the  first  two  terms  do  not  depend  on  n  but  the  third 
term  is  a  monotone  signal  with  frequency  of  2fikTs.  Figs. 
5  and  6  show  the  effect  of  chirping  rate  on  the  detection 
performance  of  the  six  integration  schemes  for  /?  =  0  Flz/s, 
37.45  Flz/s,  and  / 3  =  749  Flz/s,  respectively.  Both  the 
coherent  and  non-coherent  remain  the  same,  unchanged 
by  chirping  rate.  The  semi-coherent  for  first  lag  changed 
very  little  and  the  semi-coherent  for  first  lag  with  FFT  has 
virtually  no  change.  This  is  because  the  beat  frequency  for 
the  conjugate  products  is  2 /3kTs  =  2x749x1x0.001  =  1.5 
Flz/s  and  only  causes  a  phase  change  of  0.15  cycle  (54 
degrees)  over  the  integration  time  of  100  ms  for  k  =  1  but 
it  is  7.5  cycles  for  k  =  50.  It  therefore  has  more  effect  on 
the  semi-coherent  up  to  first  N/2  lags  as  the  purple  curve 
shifted  rightwards  in  Figs.  5  and  6. 

Flowever,  the  biggest  effect  is  on  the  practical  coherent 
with  FFT.  Because  for  the  coherent  integration  scheme, 


the  signal  frequency  changed  by  Af  =  J3NTS  = 
749x100x0.001  =  74.9  Flz  and  the  signal  energy  was 
spread  over  74.9  Flz/ 10  Flz  «  8  frequency  bins. 

Figs.  7  and  8  show  the  loss  of  integrated  signal  peak  as  a 
function  of  the  chirping  rate  for  integration  time  of  N  =  20 
and  40,  respectively.  The  three  schemes  shown,  namely, 
the  semi-coherent  for  the  first  lag  (the  blue  curve),  the 
practical  coherent  with  FFT  (the  green  curve),  and  the 
semi-coherent  up  to  first  N/2  lags  (the  red  curve),  all 
follow  the  same  sine-function  in  their  loss.  The  loss  of  the 
semi-coherent  up  to  first  N/2  lags  is  faster  than  that  of  the 
semi-coherent  for  the  first  lag  because  the  latter  only 
contains  one  beat  frequency  term  at  2 f:lTs  whereas  the 
former  contains  many  terms  at  2pkTs  for  k  >  1.  The  first 
few  nulls  of  the  .sine-function  show  up  for  the  practical 
coherent  with  FFT. 

Fig.  9  shows  the  loss  of  integrated  signal  peak  for 
integration  time  of  N  =  20,  also  as  a  function  of  the 
chirping  rate  but  expressed  in  terms  of  a  sweep  factor. 
The  sweep  factor  k  is  the  number  of  frequency  bins 
crossed  over  by  the  signal  during  the  integration  time, 
which  is  related  to  the  chirping  rate  by: 

B  =  — - — » /c=0, 1, ...,  (86) 

( TSN )2  2 

In  Fig.  9,  the  frequency  bin  is  Af  =  50  Flz  for  N  =  20.  A 
unit  sweep  factor  corresponds  to  a  chirping  rate  of  2500 
Flz/s.  Because  of  this  large  step  size,  the  nulls  are  not 
obvious  in  this  figure  as  evident  in  Fig.  7,  which  is  the 
initial  portion  of  Fig.  9.  Flowever,  the  envelope  of  signal 
peak  loss  is  clear  for  the  three  schemes. 

In  the  third  simulation,  the  input  SNR  that  is  required  to 
obtain  the  detection  performance  of  (Pd  =  0.9,  Pfa  =  0.001) 
is  determined  by  Monte  Carlo  runs  when  the  chirping  rate 
is  varied.  For  an  integration  time  of  N  samples,  the 
frequency  bin  is  l/TsN. 

Figs.  10  and  11  show  the  input  SNR  vs.  the  sweep  factor 
for  N  =  20  and  100,  respectively.  In  Fig.  10,  the  frequency 
bin  is  Af=  50  Flz.  A  unit  sweep  factor  corresponds  to  a 
chiiping  rate  of  2500  Flz/s.  As  expected  the  coherent  (the 
blue  curve),  the  non-coherent  (the  red  curve),  and  the 
semi-coherent  for  first  lag  with  FFT  (the  black  curve) 
remain  unchanged  for  different  //.  Flowever,  the  practical 
coherent  with  FFT  (the  green  curve)  crosses  (needs  a 
larger  input  SNR,  thus  having  a  worse  detection 
performance)  the  semi-coherent  up  to  first  N/2  lags  (the 
purple  curve)  at  1.8  (3750  Flz/s).  The  practical  coherent 
with  FFT  crosses  the  non-coherent  at  3  (7500  Flz/s),  the 
semi-coherent  up  to  first  N/2  lags  crosses  the  non¬ 
coherent  at  4  (10000  Flz/s),  and  the  semi-coherent  for  first 
lag  crosses  the  non-coherent  at  6  (15000  Flz/s). 

In  Fig.  11,  the  frequency  bin  is  Af=  10  Flz.  A  unit  sweep 
factor  corresponds  to  a  chirping  rate  of  100  Flz/s.  Again, 
the  coherent  (the  blue  curve),  the  non-coherent  (the  red 
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curve),  and  the  semi-coherent  for  first  lag  with  FFT  (the 
black  curve)  remain  unchanged  for  different  //.  Flowever, 
the  practical  coherent  with  FFT  (the  green  curve)  crosses 
the  semi-coherent  up  to  first  N/2  lags  (the  purple  curve)  at 
0.8  (83  Flz/s).  The  practical  coherent  with  FFT  crosses  the 
non-coherent  at  6.8  (680  Flz/s,  which  is  consistent  with 
Fig.  6),  the  semi-coherent  up  to  first  N/2  lags  crosses  the 
non-coherent  at  17  (1700  Flz/s),  and  the  semi-coherent  for 
first  lag  crosses  the  non-coherent  at  37  (3700  Flz/s).  It  is 
shown  that  the  semi-coherent  up  to  first  N/2  lags  is 
slightly  more  robust  than  the  practical  coherent  with  FFT 
against  frequency  change  as  was  previously  observed  in 
[Wirth,  2001].  At  a  first  glance,  the  semi-coherent  for  first 
lag  with  FFT  behaves  quite  similar  to  the  non-coherent  in 
both  detection  performance  and  robustness  against 
frequency  changes.  Flowever,  the  FFT  operation  provides 
an  estimate  of  the  underlying  chirping  rate,  which  is 
valuable  for  some  applications. 

Comparing  Fig.  11  to  Fig.  10  shows  that  except  for  the 
coherent  and  non-coherent  as  well  as  the  semi-coherent 


for  first  lag  with  FFT,  all  other  integration  schemes 
become  more  vulnerable  (i.e.,  experience  larger 
integration  loss)  to  changes  in  frequency  as  the  integration 
interval  gets  longer.  Yet,  longer  integration  interval 
provides  better  processing  gain  (i.e.,  integration 
improvement)  as  shown  in  Fig.  4  as  compared  to  Figs.  2 
and  3.  Flowever,  a  signal  is  more  likely  to  be  subject  to 
changes  in  frequency,  either  due  to  Doppler  frequency 
shift  or  to  clock  frequency  drift,  over  longer  intervals.  It  is 
therefore  desired  to  develop  such  an  integration  scheme 
that  is  robust  against  frequency  changes  while  providing 
larger  integration  improvements  over  long  integration 
intervals. 

One  possible  approach  is  to  extend  the  semi-coherent  for 
first  lag  with  FFT  to  other  lags.  In  other  words,  instead  of 
making  straight  summations  as  in  the  semi-coherent  up  to 
first  N/2  lags,  FFT  is  used  to  coherently  add  up  the 
conjugate  products.  This  and  other  approaches  are  under 
study  and  will  be  reported  in  future  papers  [Yang  et  al., 
2007;  2008], 


f,  =  200  Hz,  N  =  40,  t  =  1  ms 
d  s 


chirping  rate,  Hz/sec 

Fig.  8  -  Loss  of  Integrated  Signal  Peak  iovN=  40 


f .  =  200  Hz,  N  =  20,  t  =  1  ms 
d  ’  s 


f .  =  -250  Hz,  N  =  20,  runs  =  10000,  P  .  =  0.9,  Pf  =  0.001 
d  ’  ’  d  ’  fa 
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Fig.  10  -  Input  SNR  vs.  Sweep  Factor  for  N=  20 
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CONCLUSIONS 

In  this  paper,  the  coherent,  non-coherent,  and  semi- 
coherent  integration  schemes  and  their  variants  were 
presented  as  the  most  appropriate  approximation  to  the 
sufficient  statistic  of  a  log  likelihood  ratio  test  when  the 
signal  model  was  partially  known  to  different  degrees.  It 
summarized  the  formulas  for  conditional  probability 
density  functions,  probabilities  of  detection  and  false 
alarm  as  well  as  integration  improvement  factors  and 
integration  losses  for  these  integration  schemes.  Results 
of  Monte  Carlo  simulation  were  presented  to  compare  the 
detection  performance  of  these  integration  schemes  in 
terms  of  Pj  vs.  SNR  at  a  fixed  Pfa  for  signals  with  constant 
Doppler  frequency  and  for  chirp  signals.  It  also  showed 
the  curves  of  input  SNR  required  for  different  chiip  rates 
so  as  to  maintain  the  same  detection  performance.  It  can 
be  used  as  an  easy  graphic  means  to  rank  the  integration 
schemes  for  different  operation  conditions. 
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